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Abstract 

A method for determining the generahsed scahng function(s) arising in the high 
spin behaviour of long operator anomalous dimensions in the planar sl{2) sec- 
tor of = 4 SYM is proposed. The all-order perturbative expansion around 
the strong coupling is detailed for the prototypical third and fourth scaling func- 
tions, showing the emergence of the 0(6) Non-Linear Sigma Model mass-gap from 
different SYM 'mass' functions. Remarkably, only the fourth one gains contri- 
bution from the non-BES reducible densities and also shows up, as first, NLSM 
interaction and specific model dependence. Finally, the computation of the n-th 
generalised function is sketched and might be easily finalised for checks versus 
the computations in the sigma model or the complete string theory. 
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1 Introduction 



The planar s/(2) sector of A/" = 4 SYM contains local composite operators of the form 

Tr{V'Z^) + (1.1) 

where V is the (symmetrised, traceless) covariant derivative acting in all possible ways 
on the L bosonic fields Z. The spin of these operators is s and L is the so-called 'twist'. 
Moreover, this sector would be described - thanks to the AdS/CFT correspondence 
[1] - by string states on the AdSs x spacetime with AdSs and charges s and L, 
respectively. In addition, as far as the one loop is concerned, the Bethe Ansatz problem 
is equivalent to that of twist operators in QCD [21 Isj^ . 

Proper superpositions of operators (11.11) have definite anomalous dimension A de- 
pending on L, s and the 't Hooft coupling A = Sir'^g'^: 

A = L + s + 7((7,s,L), (1.2) 

where '-/{g, s, L) is the anomalous part. A great boost in the evaluation of 'y{g, s, L) in 
another sector has come from the discovery of integrability for the purely bosonic so{6) 
operators at one loop |1] . Later on, this fact has been extended to the whole theory and 
at all loops in the sense that, for instance, any operator of the form (11. ip is associated 
to one solution of some Bethe Ansatz-like equations and thus any anomalous dimension 
becomes a function, basically the energy, of one solution O [6l [7] . Nevertheless, along 
with this host of new results an important limitation emerged as a by-product of the 
on-shell (S'-matrix) Bethe Ansatz: as soon as the interaction reaches a range greater 
than the chain length, then it becomes modified by unpredicted wrapping effects. More 
precisely, the anomalous dimension is in general correct up to the L — 1 loop in the 
convergent, perturbative expansion, i.e. up to the order g^^~^ . Which in particular 
implies, - fortunately for us, - that the asymptotic Bethe Ansatz should give the right 
result whenever the subsequent limit (II. 3p is applied, as in the expansion (II. 4p . 

In fact, an important large twist and high spin scaling may be considered on both 
sides of the correspondenc^: 

s^oo, L^oo, j = - — = fixed . (1.3) 

The relevance of this logarithmic scaling for the anomalous dimension of long operators 
has been pointed out in [9j within the (one-loop) SYM theory and then in [8] and jlOj 



deep reason for that may be that one loop QCD still shows up conformal invariance, albeit 
integrability does not seem to hold in full generality (for instance, it apparently imposes aligned 
partonic helicities). 

^Actually, in string theory (semi-classical) calculations the A +oo limit needs consideration 
before any other (cf. for instance |B] and references therein), thus implying a different limit order with 
respect to our gauge theory approach (cf. below for more details). 
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within the string theory (strong couphng g ^ 1). In fact, by describing the anomalous 
dimension through a non-hnear integral equation [12] (like in other integrable theories 
[TT]). it has been recently confirmed the Sudakov leading behaviour for s — ^ +00 [TOllS] 

'j{g,s,L) = f{g,j)lns + ..., (1.4) 

thus indeed generalising to all loops a result by [9]. Actually, this statement was 
argued in [I2] by computing iteratively the solution of some integral equations and 
then, thereof, the generalised scaling function, f{g,j) at the first orders in j and g"^: 
more precisely, the first orders in g"^ have been computed for the first generalised scaling 
functions fn{g), forming 

00 

f{9,j) = Y.f-i9)r- (1-5) 

n=0 

As a by-product, the reasonable conjecture has been put forward that the two-variable 
function f{g,j) should be analytic (in g for fixed j and in j for fixed g). In [13] similar 
results have been derived for what concerns the contribution beyond the leading scaling 
function f{g) = fo{g), but with a modification which has allowed to neglect the non- 
linearity for finite L and to end-up with one linear integral equation. The latter does not 
differ from the BES one (which cover the case j = 0, cf. [7]), but for the inhomogeneous 
term, which is the sum of an integral on the one loop root density, a hole contribution 
and a known function. In this respect, we have reckoned interesting the analysis of the 
next-to-leading-order (nlo) term - still coming, for finite L, from an asymptotic Bethe 
Ansatz -, as the leading order f{g) has been conjectured to be independent of L or 
universal [S], [0] • 

For instance, this nlo term enjoys many attractive features, like its form which 
contains up to the L-linear term, i.e. fi{g)L + c{g). In particular, the linearity in L 
and the behaviour of the next terms with increasing powers (of L) has furnished us 
the inspiration for the past [H] and present calculations, respectively. In this respect, 
we deem useful to briefly introduce in the next section a suitable modification of the 
aforementioned method such that it applies to the regime fll.3p (still for any g and j). 
In fact, a suitable modification of this LIE has been already exploited and explored in 
|14] to derive still a LIE in the scaling (11.31) (for any g and j). Along this path, we want 
to determine here the generalised scaling function f{g,j) and also its constituents fn{g) 
for all the values of j and g, thus interpolating from weak to strong coupling: in fact, 
we will see in the following why the expansion around j = is suitable, efficient and 
manageable. Specifically, we will compute /«((?) for n = 2,3, 4, show that f2{g) = and 
derive the strong coupling behaviour of fs^g) and f4,{g). In this respect, a crucial point is 
the appearance in the calculation regarding f4,{g) (and onwards for n > 4) of the higher 
root densities, along with the BES density, which are not expressible in terms of the 
BES one. This novel feature (i.e. essentially new (recursive) linear integral equations 
(concerning the higher densities) beside the BES one) plays an important role in the 
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comparison with the 0(6) Non-Linear Sigma Model (NLSM). For, at the leading large g 
order, we will find perfect agreement with the infrared (i.e. strong coupling) expansion 
of the 0(6) Non-Linear Sigma Model (NLSM) [1^, thus enforcing the ultraviolet (i.e. 
weak coupling) checks and predictions by Alday and Maldacena ^U\. The IR results, 
on the one hand, explore the other regime of the asymptotically free NLSM, which 
develops a mass gap, on the other yield a convergent series (for f{g,j), in this limit, 
as a function of j and \/g in the coefficients fn{g)) [IS]. This feature, the convergence, 
likely extends to the full SYM theory, so bearing new interest to this expansion. We 
would also remark here, but we will discuss more deeply this point afterwards, that 
the agreement for f^ig) (and onwards for n > 4) is highly non-trivial since it involves 
the specific interacting theory of the 0(6) NLSM. Finally, the ideology of our method 
for computing the n-th generalised scaling function fn{g) is exposed and shown to be 
easily implementable by using analytical and numerical techniques. 

Incidentally, in the subsequent section we will also briefiy understand the corrections 
(represented by dots) to the leading Sudakov factor (11. 4p . as they are indeed formed by 
the terms like c{g) in the aforementioned nlo term. 

Note added: An interesting paper [16] appears today in the web archives. It 
seems to have some goals and equations similar to ours, although coming more directly 
from the approach [12], and giving the leading strong coupling behaviour of fsig) in 
agreement with us and confirming f2{g) = as in [12], and fi{g) as in [I4j. Nothing 
explicit we can find in it about f4,{g) and a systematic plan for all the other fn{g)- 

2 Computing the generalised scaling functions 

In the framework of integrability in A/" = 4 SYM, we have been finding useful [17] to re- 
write the Bethe equations as a unique non-linear integral equation In particular, 
we have developed (starting in pjj) a new technique in order to cope with the frequent 
(in SYM theory) case where the Bethe roots are confined within a finite interval [T3] . 
The non-linear integral equation regarding the s/(2) sector involves two functions F{u) 
and G{u,v), both satisfying linear equations. It is convenient to split F{u) into its 
one-loop and higher than one loop contributions Fq{u) and F^{u) respectively, and 
define the quantities cro(u) = -^Fq{u) and aniu) = '^F^{u), which, in the high spin 
limit, have the meaning of one loop and higher than one loop density of Bethe roots 
and holes, respectively. 

It is important to stress that in this paper we will follow the route initiated in 
[T3] . and further developed in [Tl] where the (strong coupling) behaviour of the first 
constituent fi{g) was studied in full detail. In this perspective, one of the main aims 
of the present paper is to show how we can study the generalised scaling function 

"^Only in case of nested Bethe Ansatz we may have more than one equation, but always in small 
number. 
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at all orders in j. In other words, we will show how to write down recursive linear 
integral equations for infinite many densities, cr^\ n = 0, 1, 2, . . . (the first one being 
the famous BES density [7]), each one yielding fn{g), in a way which systematically 
furnishes the n-th. one, provided those beforcl are known. Indeed, this is the main 
consequence of a suitable manipulation of the initial integral equation for the full density 
to end up with a recursive system of Fredholm (II type) linear integral equations, all 
with the same kernel (the BES one), but different inhomogeneous terms. The latter 
are the only parts entered by the previous densities and are thus responsible for the 
recursive derivation and solution. Eventually, we shall remark that this structure has 
already appeared in [HJ, allowing us to study fi{g): in fact, we did not even need 
the knowledge of the BES density a^^\ thanks to the general rule that only cr^'*, with 
m < n — 3, enter the inhomogeneous term regarding cr^'*. 

Let us now analyse in detail the linear integral equations which are satisfied by 
the one and higher loop densities and in terms of which we may write the anomalous 
dimension. 

For the time being, we want to extend to the scaling situation ( II. 3p the procedure by 
[m mj for computing the one loop observables (i.e. the physical quantities depending 
on the one loop density (To{v)). The integration on the finite range of roots, {—bo, bo), 
may be computed via an infinite integratior|f| 

bo roo 

dvf{v)ao{v) = / dvf{v)a^,{v) + 0((lns)-°°) , (2.1) 



-bo J —oo 

where the Fourier transform of the effective density cTqIv) obeys the integral equation 

r l*^! i-c ifci 

f - e 2 cos-% e~V r"^ 
<T^(fc) = -47r ^ ^ - — / rf^e^'^-XcoMcr^M -4vr5(A:)ln2, (2.2) 

2 smh y 2 smh 'y 

with (—Co, Co) the interval spanned by the internal holes [13]. This simply means 

/oo 
duXcoiu)cr'oiu) = -27r(L - 2) + 0((lns)-°°) , (2.3) 
-oo 

with the interval function Xco(^) equal to 1 if — Co < u < Cq (here the internal holes 
concentrate), and otherwise (no internal holes outside). Within the double limit (11.31) . 
the above remainders are 0((lns)-°°)@ and are exactly given by the non-linear integrals 
in [13] we have here dropped out. 

Now, we ought to briefiy comment on the corrections of the leading Sudakov for- 
mula (11.41) . namely the dots in there. We have deduced in [13] a linear integral equation 



^Starting, actually, from n — 3 backwards. 

^The parameter bo depends on the total number of roots s through the integral on (To{v) (the 
normalisation condition) which counts them (cf. [13|). 

^0{{lii s)^°°) means a remainder which goes to zero faster that any inverse power of Ins: 
lim (lns)'=0((ln.s)-°°) = 0,V/fc > 0. 
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suitable to detail the appearance of a constant term 0((lns)''), denoted in the introduc- 
tion as fi{g)L + c((?)0. Concerning this, the first addendum is the generalised scaling 
function at first j-order [H], whilst the second one gives rise to a 0((lns)°) correction 
to the leading Sudakov scaling, i.e. the j° order of some f^^\g,j) in the expansion 

^{g,s,L) = f{g,j)lns + f^'\g,j) + ... . (2.4) 

Now, similarly we may imagine that the dots should be inverse integer powers of In s, 
with coefficients, at each power, depending on g aud JE Very i„ter;sti„gly. these 
power-like corrections to the leading (11.41) seem to come out truly from our systematic 
expansion of the one-loop density, - satisfying a linear integral equation -, and its 
subsequent insertion into the inhomogeneous term of the higher loop linear integral 
equation (i.e. (12. 9p ). In fact, the last 5-term in (12.21) contributes, for instance, to 
f'^'^^g,])-, namely to the mentioned c{g), whilst it does not to /((?, j)u- In other words, 
all these terms seem to be controlled by linear equation^, after neglecting the non- 
linear integrals. Eventually, there might also be a possible non-analytic (in the variable 
Ins) correction, 0((lns)~°°). Yet, the 0((lns)~°°) term would face two problems: on 
the one side it is determined by cumbersome non-linear integrals which deny the density 
(linear) treatment, on the other side it may well become affected by the non- asymptotic 
phenomenon of wrapping. 

Therefore, as we will be constraining ourselves to the leading Sudakov factor f{g,j), 
we can neglect the aforementioned 5-term and we may, for convenience' sake, Fourier- 
transform the equations ( 12. 2p and ( 12. 3p . respectively: 

/°° dk , , , , sin kcn , , , , 

-am^^ = -2n{L-2). (2.6) 

These final relations characterise the one-loop theory in the high spin and large twist 
limit (II. 3p and need to be solved together in order to give (3"o(A;) and Cq. In specific, we 
may solve the first one, a Fredholm (II type) integral equation, by Neumann-Liouville 
(recursive) series and then expand usefully co(j) according to 



oo 



Co(j) = E^i"^^'"- (2.7) 



n=l 



^At small g, the same term has been derived in [12] by evaluating at leading order the non- linear 
integrals of a traditional NLIE approach [TT] . 

^Naturally the subsequent expansion (|2.7p and its higher loop analogue ought to be amended 
accordingly. 

^The other part of c{g) derives from the rest of the equation, i.e. the next equation (|2.5p . but it 
will be neglected in the following. 

^''This is similarly true for the linear integral formula of the anomalous dimension. 
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Already the first step of tliis procedure (tlie approximation by tlie inliomogeneous term) 
yields the first two coefficients of this series 



.(1) 



7r 

- , c, 



(2) 



-^ln2. 

4 



(2i 



Regarding the high spin behaviour of long operators, the approach to the higher 
than one loop density, crH{u), moves along similar lines. Actually, the knowledge of 
o"o(u) concurs to find it as the solution of the linear integral equation (discarding the 
following 0((lns)-°°)) 



(Jh[u 



■T ^ 1 

-iL— m 
du 



2x-{uY 



1 + 



2x+{uY 



+ - 



■ d „, . . 1 

+ i—e{u,v) + i—— 

du 1 + [u — vY 



dvxc{v) 

) 



— In 

-du 



2x+(u)x~ (v) 
2x~ {u)x+ (v) 



+ 



oo 
oo 



dv— 
du 



1 - 



dv 



In 

1 



2x+(u)x~ (v) 



1 - 



+ i9{u, v) 



2x-(u)x+(v) 



(2.9) 



vr 1 + (m — vy 
1 



°" dv 



crHiv)+ / —Xcoiv 



^ + (n-.)^ "°("^ - 



, d 
dv— 

^ du 



In 



2x+{u)x^ (v) 



+ iOiu, v) 



aH{v) + 0{{\ns) 



2x~ {u)x+ {v) 

constrained by the two conditions on c and Cq respectively 

duXc{u)[(yQ{u) + (yH{u)] = -2n{L 



duXcQ{u)aQ{u) = -27r(L - 2) 



(2.10) 



with the former fixing the interval (— c, c) where the all- loop internal holes concentrate 
as well as the latter, equivalent to (12.61) . determines the range of the one- loop internal 
holes 

As in the one loop case, in order to gain a better insight it is convenient to re-write 



^^The function 9{u, v) appearing in (|2.9p is the 'dressing factor', introduced on the string side in the 
second of [5] and finalised as meromorphic function of the couphng constant in [7]. 
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(12. 9p in terms of Fourier transforms (upon neglecting 0((lns) °^)) 



aH{k) = ttL h 



sinh ^ 



+ 



sinh 



|fc| 



dh 



-iy+^Jr{V2gk)Jr{V2gh) 



1 — sg\i{kh) _\h\_ 



e 2 + 



r=l 



+ sgTi{h)^^Cr,r+i+2u{-iy^''e ^"^^^ {^Jr-i{V2gk)Jr+2u{V2gh) 

r=2 u=0 

- Jr-i{V2gh)Jr+2u{V2gk))] \a'o{h) + anih) - 



dp 



(2.11) 

sm(h — p)c' 



_1M 
e 2 



sinh ^ i-oo 2vr 



dp . . , sin(fc - p)c e 



k-p sinh ^ i-oo 27r 



dp ^ sin(fc-p)co 



k — p 



as well as the previous conditions on the internal holes distributions 



°° dk ^„ sin kco 



dk y ^ , ^ , ...sin/uC 
— {al{k) + aHm- 



27r' 



k 



-2n{L - 2) 



-27r(L - 2) 



(2.12) 



Eventually, the three equations (12. lip and ( 12.12^ form the starting point of our all 
loop analysis, generalising the one loop statement. In fact, we may write the anomalous 
dimension (in the usual limit (11.30 ) as 



-9 



/•bo 


dv 


i 




i 


J-bo 


2^ 


_x+{ 


v) 


x^{v 


+00 


dv 






i 


oo 


2^^ 


c{v) 






+ 00 


dv 


i 




i 


oo 


2^ 


X+(f ) 


x^iy) 



[(yo{v) + aH{v)]- (2.13) 



aH{v) + 0{{\ns) 



and realise, upon comparing this with (12. lip , that 



'y{g,s,L) = -hmanik) + 0((lns) 

71 k^O 



(2.14) 



This relation extends the validity of its Kotikov-Lipatov analogue [L9\ (valid at the 
leading order Ins) to all finite orders (Ins)"", n G N. A statement of this type has been 
already noted by [I2| in connection with f{g,j) only and their equations, which have 
a different form with respect to our (12. lip , though. 



^^For instance, we verified this for f^^\g,j), albeit we do not show it here. 
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Now, a very crucial point of our construction enters the stage. The key equation 
(12. lip concerning the all loop density aniu) and also the conditions (12.121) show up a 
deep, even conceptual, difficulty to be solved for generic values of the parameters g and 
j (or s, L). Yet, finding the solution becomes easier if we think of it and of the internal 
holes boundary c(j) as an expansion in (non- negative) powers of j 

oo oo 

aniu) = [J2 crP{u)r] In s , c(j) = ^ c(")j- . (2.15) 

n=0 n=l 

For in this way the equation (12. lip breaks down into a recursive linear system of integral 
equations for the densities cr^'', each one yielding /„((?) according to (I2.14p : 

Mg) = -^m/-\k). (2.16) 

Importantly, each o"^^ is governed by a Fredholm (II type) linear integral equation with 
always the same kernel (indeed the BES one [7]), but a different inhomogeneous term. 
This contains a part equal for any n and given by the one loop density (and an additional 
known function) and another part which instead varies, but is still expressible - after 
using the conditions ( 12.12P - in terms of only (t^^\ with m < n—3. Therefore, a recursive 
solution procedure can be contrived, at least in principle. In the ffist step of the 
recursive structure has already appeared, without any reference to the BES density 
afp , and allowed us to analyse fi{g)- From a physical point of view, the expansion in 
j means that we are focusing our attention on the regime j <^ g, with fixed coupling g. 
Unfortunately, semi-classical string results concern a different regime where g is very 
large and then j scales, instead, accordingly, i.e. j ^ g (in fact, the scaling variable 
y = j/g is fixed) [8]. Nevertheless, our regime is suitable for comparison with the string 
reduction to the 0(6) NLSM which happens for j <^ g as brilliantly conjectured by 
Alday and Maldacena on a geometric ground [lOj. These authors have also checked 
their statement by using the perturbative renormalisation procedure in the UV regime, 
where the mass-gap m{g) is much smaller than the energy density j, i.e. j ^ fnig)- On 
the contrary, our kind of expansion obliges us to consider j small and then to consider 
the other NLSM regime, the IR one j <^ fn{g)- In this case, the perturbative field 
theory methods are much less effective, but fortunately the Bethe Ansatz formulation, 
based on the S'-matrix, furnishes easily the necessary data still in a systematic way [TJ] . 
The first check of the mass gap arising has been confirmed for the first time in [13] by 
computing fi{g), while it could not be shown exactly in the strong coupling behaviour 
of the cusp anomalous dimension fo{g) = f{g) because of the UV cut-off effects on the 
NLSM vacuum energy. Instead, f\{g) gives rise to the simplest manifestation of the 
mass-gap and the other scaling functions fn{g) - as we are going to see - present other 
different functions (of g), which all coincide with the mass-gap in the strong coupling 
regime. 
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In specific, we will illustrate explicitly our general construction in the following 
section giving as prototypical example the first fully general case, i.e. that concerning 
fi{g). As a matter of fact, this is the lowest order at which this computation is not 
reducible to the knowledge of the BES density o"^'* (at least by us and by now), as it 
shows up the explicit appearance of the higher order density <J^^ ^ in the inhomogeneous 
term. Curiously, it is also the lowest order at which the comparison with the 0(6) 
NLSM becomes constraining, since right at this order the specific interaction starts 
appearing. In other words, up to the order j^, the strong coupling expansion obtained 
for the 0(6) NLSM only relies on the non- interacting Fermi gas theory, which has no 
information about the interaction and is in some sense " universal" , being , in particular, 
the same for any of the 0{N). In other words, the exact leading order of fn{g) ^ as 
it comes out from the gauge theory - is in perfect agreement with the corresponding 
one as worked out in [ISj within the 0(6) NLSM. In the following, we have decided 
to illustrate the computation up to the first general example, fi{g), which teaches us 
the general procedure. Nevertheless, all the details and important subtleties about 
the latter need to be given extensively in a separate publication; for this reason we 
have simplified the illustration and have limited ourselves to an effective numerical 
presentation of the leading non-perturbative terms, which are responsible for the mass- 
gap. Yet, the different physical origin of the fn{g) will clearly prove that the complete 
NLSM allows (with its unique mass-gap) only for their leading terms, though giving 
rise to an impressive exact agreement. 

3 On the calculation of the generalised scaling func- 
tions: a sketch 

In the present section we will sketch how to compute the coefficients of the expansion 
( I2.15p . so that, by means of the equality fl2.14p . we can gain all the generalised scaling 
functions fn{g)- In this respect, computing the first orders in the series for c(j) fl2.15p 
is instructive to see how the procedure works and will turn out useful for expanding the 
subsequent forcing terms (cf. sub-section 3.3). The computation goes along the similar 
lines of the one loop theory and yields 

cC-.^L^, e>"^- /'°^-f<°' . (3.1) 
3.1 The second generalised scaling function 

We now show that o"^ (m) = 0, so that we obviously have f2{g) = 0. Let us consider 
the r.h.s. of fl2.9p . The first term is clearly proportional to j Ins, so it does not appear 
in the equation for a^^\u). The second and the fifth term both have the form, with two 
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different functions f{v) 



+00 



dvf{v)cr{v)xd{v) 



dp ^ , ,sm(k — p)d 



27r 



k — p 



(3.2) 



where a, d stand for a^ + an, c if we consider tlie second term and for ctq, cq if we consider 
the fifth term. Using the normalization condition 



one can show that 



dp ^ , .sinpd „ ^ ^ , 
^a{p) — — = -27rL + O Ins 

/TT p 



'^'^ dp ^ , ,sm(k — p)d r r. ■ ^/7S\Ti 

(y{p) — ) — = [-271 J + 0{d^)]lns. 



2tx 



p 



(3.3) 



(3.4) 



Since d starts from order j in its expansion, the second and the fifth term in the r.h.s. 
of (12.91) lack of the order In s terms in their expansion. The same reasoning, apphed 
to the second term in the rhs of (12.51) - the one containing the integral - implies that 
also this term lacks of the order In s. Therefore, the third term in the r.h.s. of (12. 9p 
is missing the quadratic order as well. It follows that the equation for a\j {u) is 

' (3,5) 



.(2)/ 



TT 1 + (m — f )^ 

1 - 



dv-^ 
du 



In 



2x+{u)x~ (v) 



+ i6{u, V) 



whose solution is, of course, aff(u] 



2x-(u)x+(v) . 

0. Therefore f2{g) = 0. 



a 



(2), 
H I 



3.2 The third generaUsed scahng function 

In this sub-section, we concentrate our attention on the equation for a^^\u). We find 
convenient to pass to the Fourier transforms of all the involved quantities. Let us define 
the even function 

^''\k) = ^-^^aP{k), (3.6) 



which has the property 



27r\k\ 



hms^'\k) = -Ug) 



For future convenience, we will study the new function 



.(3) 



TT 



(k)--\k\e~^ 



1 

16 



[a(o)(0)]^ 



(3.7) 



(3.^ 



which depends on the all loops density of roots at (the zero-th) order Ins (satisfying 
the BES equation [7]) at m = (real space) 

a(o)(0) = -4 + ag)(0), (3.9) 
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(being the zero-th one-loop density in zero a^\o) = —4) and still keeps the key property 
fl2.16p in the form 

hmS'^'\k) = lMg). (3.10) 
Let us also introduce the functions 

r^^^^. (3.11) 



6[-4 + ag)(0)]2 io sinhf ' 

and restrict our analysis within the domain k > 0. Upon expanding in series of Bessel 
functions (the so-called Neumann's expansion) as 

S(=)(i) = ± sg'(g)M^ + ^ SllM '^'-^'f""^ ■ (3.12) 

p=l p=l 

introducing the quantified 

z„,„(,)^rf^^i»^. (3.13) 

Jq n e — i 

we finally obtain the following infinite system of linear equations: 

oo oo 
SSig) = 4?(^7)-4p^Z2,,2„^(^?)5S(^?)+4p5^Z2,,2™-l(^?)S(l,((^^ 

m=l m=l 

oo 

sS-iia) = A^iU9)-2{2p-l)J2Z2p-iM9)S?l{9)- (3.14) 

m=l 

oo 

- 2(2p -1)Y, Z2,^i,2,n^i{9)S^il_,{g) . 



m=l 



The Neumann's decomposition fl3.12p reveals itself useful also thanks to the special role 
of the first component in the key relation 

f3{9) = V2gS['\g), (3.15) 

easily derived from (13.101) . 

As a first approach, if we forget the non-perturbative contributions around g = 
-I-OC0, we may seek for the solution to the equations (13.140 as an asymptotic series 

oo Q(3,2n) oo Q(3,2n-1) 

n=l ^ n=2 ^ 



^•^Their appearance comes out, also but not only, from the coefficients of the dressing phase Cr^sid) = 
2cos[f(s-r - 1)] {r-l){s~l)Zr-i,s-i{g) as in [7]. 

^"^Sometimes these are improperly called non- analytic, in the sense that they have this character, 
although the following series is (as a consequence!) not convergent at all. 
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The first two coefficients of fl3.16p at tlie orders 1/g and 1/ g'^ have the same structure as 
the corresponding ones in the case of the ffist constituent fi{g) [H]. Therefore, we are 
naturally led to make the same Ansatz on the form of the coefficients entering fl3.16p : 

e{3,2n) ^ r(m + n) iNl+nT(3) ^ . ^ . 

= ^^ Y{m -n+1) ^" ' ' ' 

(3.17) 

siZV^ = i2rn-l) l^^^_^^^^^^ n>2, m>2. 

Analogously to the case of the ffist constituent, this Ansatz implies that S^^^i and 
S^2m^^ are different from zero only if n < m. The coefficients foi^^ are unknown and they 
are determined by inserting the Ansatz fl3.16[ 13.171) into the equations (13.141) . 
In this way, implementing the elegant asymptotic expansion 

we obtain that the coefficients hn satisfy the following (infinite) linear system 
62 = g(-l)-.2-4^6£U„,,, n>l, 

m=0 ^ '' 



(3.19) 

^^"+1 - 6[a(0)(0)]2^^"^^ (2n)! ^2n+2^(-l) 2 ^^&2n-2..+2, n>l. 



m=0 



Again, the solution to these equations comes from the comparison with the correspond- 
ing ones regarding the ffist generalised scaling function [H], in the form of the simple 
mapping 



00 

b^'\t) = J2b^^k\ (3.21) 



Still the generating function 



is simple, though, namely 



^^'^^^^ " 6[a(o)(0)]2 cos - sin ^ ■ ^^'^^^ 

This concludes what concerns the asymptotic part of the study about fsi^g). Instead, 
the non-perturbative (or non-analytic) terms will be studied in subsection 3.4 (in com- 
parison with those for /4 ((?)), as they give naturally rise to the limiting mass-gap. 



n=2 



2 2f 
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3.3 On the fourth and the other generahsed scahng functions 

Thanks to the lesson of the previous sub-section, we are now in the position to formulate 
a general scheme for computing the n-th generalised scaling function fn{g) for n > 2 
at arbitrary value of the coupling constant. In particular, we can show how the first 
computation, which is not reducible to the leading BES equation, is exactly that for 
fi{g). Nevertheless, we will show how the computation for fi{g) may be reduced to 
that for fz{g)-, via a proportionality factor. The latter turns out to be the only part 
which is not reducible to the BES equation. Actually, similar situations would appear 
for the other scaling functions, as may be guessed from the results of this sub-section. 

In the hmit (11. 3p the scaling may be thought of as governed by the function S{k), 
defined through 



(3.23) 

and with the crucial expansion 

00 

S{k) = J2S^''\k)r- (3-24) 

n=0 

Let us now concentrate our attention on S^"\k), with n > 2, since the case n = 1 is 
slightly different from the generic one (as shown in [H]). Restricting the domain to 
> 0, we expand again in Neumann's series (of Bessel functions) 

S^-\k) = ± S^;\g)^^^^^y^ + f: SiUg /''~'^f''^ , (3.25) 
p=i p=i 

and find the following equations for the coefficients of S^'^\k), with n >2, 

00 00 
Si;\g) = 4?(^?)-4p$^Z2,,2™(<7)5S(^?)+4p^Z2,,2™_i(^)41,(^) 

m=l m=l 

00 

Si;Ug) = 4?-i(^7)-2(2p-l)5^Z2,-i,2™(^7)4:](^)- (3.26) 

m=l 

00 

- 2{2p-i)J2z2,-i,2m-i{g)st2-i{g)- 



m=l 

The forcing terms Ai^'' are given by the following integrals 



(3.27) 



where the symbol |jn means that we wish to keep only the coefficient of and also 
neglect its overall factor Ins. 
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Thanks to fl2.16p and also by making use of (12.121) . we can generalise (13.101) to 

limS^-\k) = lUg). (3.28) 

As before (cf. (I3.15P ). only the first component enters the expression for the generalised 
constituent 

Uig) = V2gS[''\g). (3.29) 

From the relations (I2.121I3TT]) . we can gain the relevant power series expansion (here 
limited at the order for simplicity's sake) 



oo 



oo 



|^2!l^^^^K(p) + *„(p)l= (3.30) 
Zn ri — p 



3|-4 + ag'(0)|= 3 |-4 + a2'(0)| 



3 



In s . 



From direct inspection of the second and third term in the r.h.s. of (13.301) . we easily 
realise the proportionality between f^i^g) and f3{g): 

which allows us to use the results about /s, provided the non-BES pre-factor be com- 
puted. In fact, this pre-factor contains explicitly the first correction to the BES ap- 
proach, i.e. a^^\o), whose computation cannot apparently be derived from the BES 
equation. 

In general, it is thus possible to systematically expand the solution (I3.23P into the 
form (13. 24113.261) . namely order by order in j up to the desired order. Albeit this 
procedure seems to be in principle straightforward, it contains many intrigued and 
intriguing details which are to be left for a dedicated publication. Therefore, we plan to 
analyse in the next section the non-analytic contributions in a numerical fashion (as we 
reckon it quite effective for the current presentation), leaving the analytic study about 
them for that publication. 

4 The non-analytic contributions from SYM and 
the 0(6) NLSM from string dual 

Let us show how it is possible to make contact with the results of the 0(6) NLSM 
[T3] , which appears in the dual string theory description, by analyzing the non-analytic 
exponential contributions to fn{g)- 



^^The term proportional to j in (|3.30p needs careful consideration since it receives additional con- 
tributions to produce the final inhomogeneous equation of [f4] . 
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First of all, let us focus our attention on f^{g). From the form of the inhomogeneous 
term fl3.1ip (given by the coefficient of (13.301) ). we may think of introducing a reduced 
form of the equations (13.141) 

oo oo 

S'S^''\9) = AfJ'-\g) - 4p 5^ Z,,M9)SiT'i9) + ^pJ2 ^2p,2™-i(^7)4S':f (^7) 

m=l m=l 

oo 

m=l 

oo 

- 2{2p -1)Y, Z2,-i,2n.-i{9)S^il':f{9) , 



m=l 

in=l 

by defining a new inhomogeneous term 



n)-'-^^^^f^m.)-r r iKu'-^, (4,2) 

TT^ Jq smh " 



which, thanks to linearity, easily entails 



7r2 



Actually, the reduced system may be introduced for any n of the systems (13.261) and 
their peculiarity will be explained in a future publication, as they admit a little involved 
mapping into the BES equation (only as far as the computation of the first component or 
fn{9) is concerned). On the contrary, the mapping could not be found for the complete 
systems (I3.26P (but, of course, for n = 1,2,3). Thanks to the proportionality (13.311) . 
the reduced part of f^ig) is still a very simple matter, as it virtually coincides with 
fsig), but the prefactor with (Jh\o) could not be reformulated in terms of the BES 
equation: this makes f^ig) the first (simple) general case. 

Besides its analytical extra-value, the matrix inhomogeneous (linear) system ( 14. ip 
- indeed equivalent to the linear integral equation at this order - may be explored 
numerically at first quick instance, by truncating the dimension of the matrix kernel 
(and vectors). Naturally, this treatment turns out into the lines by [20] for the BES 
case, as it shares the same (matrix) kernel: this strategy has been already efficiently 
exploited in [Hj for the computation of the first mass-gap, namely that related to fi{g)- 
Thus, we obtain a new mass-gap behaviour with the leading form 

fr'{g) = kl''9'^'e-V^^'. (4.4) 



A fortiori, the numerical method [20] for the BES case in the Fourier space may be 
directly applied to obtain the real space BES estimate 

a(;)(0) =4 + fcfV/V^^ (4.5) 
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along with the best fit estimates for ki and k2 

fcf* = -7.1166 ±0.0005, A;f* = 5.5896 ±0.0005. (4.6) 

For clarity's sake, we should recall the exact mass-gap formula for the 0(6) NLSM via 
a convergent Taylor series of the 't Hooft coupling around 1/g = 



m 



(^7)=A:,^/V^^fl±^±^±...) , (4.7) 



with the embedding pre-factor, fc, fixed by the weak-coupling perturbation theory (upon 
comparing versus the one- loop string result [S]) by [TD] 

25/8^1/4 

Eventually, we can easily realise how strongly the following exact values for ki, k2 
become suggested: 

fci = -vrfc, k2 = —k. (4.9) 

To conclude the part on fsig), we may write down its numerical expression at strong 
coupling 

2 

fJg) = — (0.110366 ± 0.000089) g'^^^e^'^ , (4.10) 
6 

or, alternatively, making use of the guessed exact values for ki and /c2, the exact (strong 
coupling) value 

Mg) = T^X^y (4.11) 

This is the same result as that derived by the exact expansion of the 0(6) NLSM energy 
in the strong regime j <^ rn{g) as presented in [15], though this physics is explained 
uniquely by the free theory. 

Therefore, we ought to analyse in more detail the situation regarding fiig) to gain 
some physical insight into this peculiar limit. For this purpose we shall evaluate the 
ratio fi{g) I f^{g), as coming from (13.311) . at large values of g. 

Let us first compute analytically the strong coupling limit of (0) (indeed, a non- 
vanishing constant, as we will see below). Simply using the definition of S'^^\k) as 
Neumann's expansion with coefficients S^{g) [14j, we obtain 



^^AU the other coefficients of the series, i.e. ai, 02, . . . become fixed by the string UV embedding, in 
fact realised by the massive excitations. The latter are also responsible for the quenched exponential 
terms, 0(e~^^^), allowed by the SYM theory, and instead forbidden by the NLSM. 
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As we have recently found the exact asymptotic solution (around g = +00) for the 
coefficients S"™"* (g) [H] , we only need to sum over the index m so that 



dk 



sinh I 



n=l n=l 



(4.13) 



Eventually, we know from (4.4) of [H] the generating function h^^\k) and thus we may 
integrate to obtain at leading order the constant value 



aJ!Vo) = 2 / dx— 



^ cosh X sinh x 



cosh 2x cosh 2x 



31n2-^. (4.14) 



On the other hand, we have already obtained above (cf. the discussion concerning 
fsig)) the leading behaviour 



-4 + aP{0) = ~7cm{g), (4.15) 

where mlg) is exactly the 0(6) NLSM mass gap (the dependence on the coupling g by 
the density values in zero is omitted). Therefore, we can conclude that 



Ujg) _ 21n2 + 7r 

fsig) T^Mg) 

which implies an exact relation for the first strong coupling term 



(4.16) 



^ , , 7r(21n2 + 7r) ,^ 
^-(^^ = - 2A^^ig) • ^'-''^ 

This prediction is in remarkable agreement with the computation of the ground state 
energy of the 0(6) NLSM [15] and extends the agreement between the latter model and 
the gauge theory up to the first order (j^) where the interaction goes on the stage. On 
the contrary, the 0(6) NLSM approximation would no longer be applicable for the next 
(exponential) order, 0(e ^ ), present in the SYM theory with all the subsequent 
ones, as it does not show up any terms of this form. In fact, all the exponential 
corrections should be generated by the massive (fermionic and bosonic) excitations of 
the dual string. Hence, due to their origin, these corrections should be different in the 
different fn{g)'- for the time being we can state the difference at the next-to-leading 
order among the mass terms coming from /i, and /4 in the gauge theory. A string 
theory verification of this departure from the NLSM regime would be desirable. 



5 Summary 

In this work we have initiated a general method for computing all the generalised scaling 
function fn{g) appearing in the sl{2) sector of A/" = 4 SYM when both the spin s and 

^^Not to be confused with the mass-gap m{g). 
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the twist L are very large, while the ratio j = L/lns stays fixed (namely the limit 
(11.3p ) . Our method relies on solving recursively the linear equation (12. 9p for the higher 
loop density of Bethe roots aniu), which in its turn is completely determined by the 
analogous solution of the one loop linear equation for the density o"o(m). In particular, 
we have focused our attention on the third and the fourth constituents to disentangle the 
emergence, in the SYM theory, of different 'mass' terms, which all flow to the mass-gap 
of the 0(6) Non-Linear Sigma Model: comparison with the full string theory would be 
highly interesting, though still missing. Very peculiarly, the convergence of the ground 
state energy as a Taylor series in j in the NLSM [15] strongly suggests how its extension 
to the full theory, f{g,j) = Yl'^=o fn{9)j^, should be convergent too. In conclusion, we 
have left apart some details about the analytic calculations and the general structure, 
bearing in mind a more complete and systematic study for the near future. 
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